Replication of the association of common rs9939609 variant of FTO with increased BMI in an Australian adult twin population but no evidence for gene by environment (G x E) interaction

OBJECTIVE: To further investigate a common variant (rs9939609) in the fat mass- and obesity-associated gene (FTO), which recent genome-wide association studies have shown to be associated with body mass index (BMI) and obesity. DESIGN: We examined the effect of this FTO variant on BMI in 3353 Australian adult male and female twins. RESULTS: The minor A allele of rs9939609 was associated with an increased BMI (P=0.0007). Each additional copy of the A allele was associated with a mean BMI increase of approximately 1.04 kg/m(2) (approximately 3.71 kg). Using variance components decomposition, we estimate that this single-nucleotide polymorphism accounts for approximately 3% of the genetic variance in BMI in our sample (approximately 2% of the total variance). By comparing intrapair variances of monozygotic twins of different genotypes we were able to perform a direct test of gene by environment (G x E) interaction in both sexes and gene by parity (G x P) interaction in women, but no evidence was found for either. CONCLUSIONS: In addition to supporting earlier findings that the rs9939609 variant in the FTO gene is associated with an increased BMI, our results indicate that the associated genetic effect does not interact with environment or parity.

Gene by environment QTL mapping through multiple trait analyses in blood pressure salt-sensitivity: identification of a novel QTL in rat chromosome 5

Background The genetic mechanisms underlying interindividual blood pressure variation reflect the complex interplay of both genetic and environmental variables. The current standard statistical methods for detecting genes involved in the regulation mechanisms of complex traits are based on univariate analysis. Few studies have focused on the search for and understanding of quantitative trait loci responsible for gene × environmental interactions or multiple trait analysis. Composite interval mapping has been extended to multiple traits and may be an interesting approach to such a problem. Methods We used multiple-trait analysis for quantitative trait locus mapping of loci having different effects on systolic blood pressure with NaCl exposure. Animals studied were 188 rats, the progenies of an F2 rat intercross between the hypertensive and normotensive strain, genotyped in 179 polymorphic markers across the rat genome. To accommodate the correlational structure from measurements taken in the same animals, we applied univariate and multivariate strategies for analyzing the data. Results We detected a new quantitative train locus on a region close to marker R589 in chromosome 5 of the rat genome, not previously identified through serial analysis of individual traits. In addition, we were able to justify analytically the parametric restrictions in terms of regression coefficients responsible for the gain in precision with the adopted analytical approach. Conclusion Future work should focus on fine mapping and the identification of the causative variant responsible for this quantitative trait locus signal. The multivariable strategy might be valuable in the study of genetic determinants of interindividual variation of antihypertensive drug effectiveness.

A deleterious gene-by-environment interaction imposed by calcium channel blockers in Marfan syndrome.

Calcium channel blockers (CCBs) are prescribed to patients with Marfan syndrome for prophylaxis against aortic aneurysm progression, despite limited evidence for their efficacy and safety in the disorder. Unexpectedly, Marfan mice treated with CCBs show accelerated aneurysm expansion, rupture, and premature lethality. This effect is both extracellular signal-regulated kinase (ERK1/2) dependent and angiotensin-II type 1 receptor (AT1R) dependent. We have identified protein kinase C beta (PKCβ) as a critical mediator of this pathway and demonstrate that the PKCβ inhibitor enzastaurin, and the clinically available anti-hypertensive agent hydralazine, both normalize aortic growth in Marfan mice, in association with reduced PKCβ and ERK1/2 activation. Furthermore, patients with Marfan syndrome and other forms of inherited thoracic aortic aneurysm taking CCBs display increased risk of aortic dissection and need for aortic surgery, compared to patients on other antihypertensive agents.

Molecular Genetic Basis of Opposite Gene by Environment Interactions in Reading Disability and Attention Deficit/Hyperactivity Disorder

The goal of this study is to better understand the genetic basis of Reading Disability (RD) and Attention Deficit Hyperactivity Disorder (ADHD) by examining molecular G x E interactions with parental education for each disorder. Research indicates that despite sharing genetic risk factors, RD and ADHD are influenced by different types of G x E interactions with parental education - a diathesis stress interaction in the case of ADHD and a bioecological interaction in RD. In order to resolve this apparent paradox, we conducted a preliminary study using behavioral genetic methods to test for G x E interactions in RD and the inattentive subtype of ADHD (ADHD-I) in the same sample of monozygotic and dizygotic Colorado Learning Disabilities Research Center same-sex twin pairs (DeFries et al., 1997), and our findings were consistent with the literature. We posited a genetic hypothesis for this opposite pattern of interactions, which suggests that only genes specific to each disorder enter into these opposite interactions, not the shared genes underlying their comorbidity. This study sought to further investigate this paradox using molecular genetics methods. We examined multiple candidate genes identified for RD or related language phenotypes and those identified for ADHD for G x E interactions with parental education. The specific aims of this study were as follows: 1) partition known risk alleles for RD and/or related language phenotypes and ADHD-I into those which are pleiotropic and non-pleiotropic by testing each risk allele for association with both RD and ADHD-I, 2) explore the main effects of parental education on both RD and ADHD-I, 3) address G-E correlations, and 4) conduct exploratory G x E interaction analyses in order to test the genetic hypothesis. Analyses suggested a number of pleiotropic genes that influence both RD and ADHD; however, results did not remain after correcting for multiple comparisons. Although exploratory G x E interaction findings were not significant after multiple comparison correction, results suggested a G x E interaction in the bioecological direction with KIAA0319, parental education, and ADHD-I. Given the limited power in the current study, replication of these findings with larger samples is necessary.

Gene by smoking interaction in blood pressure

Hypertension (HT) is mediated by the interaction of many genetic and environmental factors. Previous genome-wide linkage analysis studies have found many loci that show linkage to HT or blood pressure (BP) regulation, but the results were generally inconsistent. Gene by environment interaction is among the reasons that potentially explain these inconsistencies between studies. Here we investigate influences of gene by smoking (GxS) interaction on HT and BP in European American (EA), African American (AA) and Mexican American (MA) families from the GENOA study. A variance component-based method was utilized to perform genome-wide linkage analysis of systolic blood pressure (SBP), diastolic blood pressure (DBP), and HT status, as well as bivariate analysis for SBP and DBP for smokers, non-smokers, and combined groups. The most significant results were found for SBP in MA. The strongest signal was for chromosome 17q24 (LOD = 4.2), increased to (LOD = 4.7) in bivariate analysis but there was no evidence of GxS interaction at this locus (p = 0.48). Two signals were identified only in one group: on chromosome 15q26.2 (LOD = 3.37) in non-smokers and chromosome 7q21.11 (LOD = 1.4) in smokers, both of which had strong evidence for GxS interaction (p = 0.00039 and 0.009 respectively). There were also two other signals, one on chromosome 20q12 (LOD = 2.45) in smokers, which became much higher in the combined sample (LOD = 3.53), and one on chromosome 6p22.2 (LOD = 2.06) in non-smokers. Neither peak had very strong evidence for GxS interaction (p = 0.08 and 0.06 respectively). A fine mapping association study was performed using 200 SNPs in 30 genes located under the linkage signals on chromosomes 15 and 17. Under the chromosome 15 peak, the association analysis identified 6 SNPs accounting for a 7 mmHg increase in SBP in MA non-smokers. For the chromosome 17 linkage peak, the association analysis identified 3 SNPs accounting for a 6 mmHg increase in SBP in MA. However, none of these SNPs was significant after correcting for multiple testing, and accounting for them in the linkage analysis produced very small reductions in the linkage signal. ^ The linkage analysis of BP traits considering the smoking status produced very interesting signals for SBP in the MA population. The fine mapping association analysis gave some insight into the contribution of some SNPs to two of the identified signals, but since these SNPs did not remain significant after multiple testing correction and did not explain the linkage peaks, more work is needed to confirm these exploratory results and identify the culprit variations under these linkage peaks. ^

PPARGC1A sequence variation and cardiovascular risk-factor levels: a study of the main genetic effects and gene x environment interactions in children from the European Youth Heart Study.

AIMS/HYPOTHESIS: The PPARGC1A gene coactivates multiple nuclear transcription factors involved in cellular energy metabolism and vascular stasis. In the present study, we genotyped 35 tagging polymorphisms to capture all common PPARGC1A nucleotide sequence variations and tested for association with metabolic and cardiovascular traits in 2,101 Danish and Estonian boys and girls from the European Youth Heart Study, a multicentre school-based cross-sectional cohort study. METHODS: Fasting plasma glucose concentrations, anthropometric variables and blood pressure were measured. Habitual physical activity and aerobic fitness were objectively assessed using uniaxial accelerometry and a maximal aerobic exercise stress test on a bicycle ergometer, respectively. RESULTS: In adjusted models, nominally significant associations were observed for BMI (rs10018239, p = 0.039), waist circumference (rs7656250, p = 0.012; rs8192678 [Gly482Ser], p = 0.015; rs3755863, p = 0.02; rs10018239, beta = -0.01 cm per minor allele copy, p = 0.043), systolic blood pressure (rs2970869, p = 0.018) and fasting glucose concentrations (rs11724368, p = 0.045). Stronger associations were observed for aerobic fitness (rs7656250, p = 0.005; rs13117172, p = 0.008) and fasting glucose concentrations (rs7657071, p = 0.002). None remained significant after correcting for the number of statistical comparisons. We proceeded by testing for gene x physical activity interactions for the polymorphisms that showed nominal evidence of association in the main effect models. None of these tests was statistically significant. CONCLUSIONS/INTERPRETATION: Variants at PPARGC1A may influence several metabolic traits in this European paediatric cohort. However, variation at PPARGC1A is unlikely to have a major impact on cardiovascular or metabolic health in these children.

From nature versus nurture, via nature and nurture, to gene x environment interaction in mental disorders

It is now generally accepted that complex mental disorders are the results of interplay between genetic and environmental factors. This holds out the prospect that by studying G x E interplay we can explain individual variation in vulnerability and resilience to environmental hazards in the development of mental disorders. Furthermore studying G x E findings may give insights in neurobiological mechanisms of psychiatric disorder and so improve individualized treatment and potentially prevention. In this paper, we provide an overview of the state of field with regard to G x E in mental disorders. Strategies for G x E research are introduced. G x E findings from selected mental disorders with onset in childhood or adolescence are reviewed [such as depressive disorders, attention-deficit/hyperactivity disorder (ADHD), obesity, schizophrenia and substance use disorders]. Early seminal studies provided evidence for G x E in the pathogenesis of depression implicating 5-HTTLPR, and conduct problems implicating MAOA. Since then G x E effects have been seen across a wide range of mental disorders (e.g., ADHD, anxiety, schizophrenia, substance abuse disorder) implicating a wide range of measured genes and measured environments (e.g., pre-, peri- and postnatal influences of both a physical and a social nature). To date few of these G x E effects have been sufficiently replicated. Indeed meta-analyses have raised doubts about the robustness of even the most well studied findings. In future we need larger, sufficiently powered studies that include a detailed and sophisticated characterization of both phenotype and the environmental risk.

The effect of genotype-by-environment interaction on longitudinal triglyceride profiles in the Bogalusa Heart Study

Triglyceride levels are a component of plasma lipids that are thought to be an important risk factor for coronary heart disease and are influenced by genetic and environmental factors, such as single nucleotide polymorphisms (SNPs), alcohol intake, and smoking. This study used longitudinal data from the Bogalusa Heart Study, a biracial community-based survey of cardiovascular disease risk factors. A sample of 1191 individuals, 4 to 38 years of age, was measured multiple times from 1973 to 2000. The study sample consisted of 730 white and 461 African American participants. Individual growth models were developed in order to assess gene-environment interactions affecting plasma triglycerides over time. After testing for inclusion of significant covariates and interactions, final models, each accounting for the effects of a different SNP, were assessed for fit and normality. After adjustment for all other covariates and interactions, LIPC -514C/T was found to interact with age3, age2, and age and a non-significant interaction of CETP -971G/A genotype with smoking status was found (p = 0.0812). Ever-smokers had higher triglyceride levels than never smokers, but persons heterozygous at this locus, about half of both races, had higher triglyceride levels after smoking cessation compared to current smokers. Since tobacco products increase free fatty acids circulating in the bloodstream, smoking cessation programs have the potential to ultimately reduce triglyceride levels for many persons. However, due to the effect of smoking cessation on the triglyceride levels of CETP -971G/A heterozygotes, the need for smoking prevention programs is also demonstrated. Both smoking cessation and prevention programs would have a great public health impact on minimizing triglyceride levels and ultimately reducing heart disease. ^

BAYESIAN STATISTICAL METHODS IN GENE-ENVIRONMENT AND GENE-GENE INTERACTION STUDIES

Complex diseases such as cancer result from multiple genetic changes and environmental exposures. Due to the rapid development of genotyping and sequencing technologies, we are now able to more accurately assess causal effects of many genetic and environmental factors. Genome-wide association studies have been able to localize many causal genetic variants predisposing to certain diseases. However, these studies only explain a small portion of variations in the heritability of diseases. More advanced statistical models are urgently needed to identify and characterize some additional genetic and environmental factors and their interactions, which will enable us to better understand the causes of complex diseases. In the past decade, thanks to the increasing computational capabilities and novel statistical developments, Bayesian methods have been widely applied in the genetics/genomics researches and demonstrating superiority over some regular approaches in certain research areas. Gene-environment and gene-gene interaction studies are among the areas where Bayesian methods may fully exert its functionalities and advantages. This dissertation focuses on developing new Bayesian statistical methods for data analysis with complex gene-environment and gene-gene interactions, as well as extending some existing methods for gene-environment interactions to other related areas. It includes three sections: (1) Deriving the Bayesian variable selection framework for the hierarchical gene-environment and gene-gene interactions; (2) Developing the Bayesian Natural and Orthogonal Interaction (NOIA) models for gene-environment interactions; and (3) extending the applications of two Bayesian statistical methods which were developed for gene-environment interaction studies, to other related types of studies such as adaptive borrowing historical data. We propose a Bayesian hierarchical mixture model framework that allows us to investigate the genetic and environmental effects, gene by gene interactions (epistasis) and gene by environment interactions in the same model. It is well known that, in many practical situations, there exists a natural hierarchical structure between the main effects and interactions in the linear model. Here we propose a model that incorporates this hierarchical structure into the Bayesian mixture model, such that the irrelevant interaction effects can be removed more efficiently, resulting in more robust, parsimonious and powerful models. We evaluate both of the 'strong hierarchical' and 'weak hierarchical' models, which specify that both or one of the main effects between interacting factors must be present for the interactions to be included in the model. The extensive simulation results show that the proposed strong and weak hierarchical mixture models control the proportion of false positive discoveries and yield a powerful approach to identify the predisposing main effects and interactions in the studies with complex gene-environment and gene-gene interactions. We also compare these two models with the 'independent' model that does not impose this hierarchical constraint and observe their superior performances in most of the considered situations. The proposed models are implemented in the real data analysis of gene and environment interactions in the cases of lung cancer and cutaneous melanoma case-control studies. The Bayesian statistical models enjoy the properties of being allowed to incorporate useful prior information in the modeling process. Moreover, the Bayesian mixture model outperforms the multivariate logistic model in terms of the performances on the parameter estimation and variable selection in most cases. Our proposed models hold the hierarchical constraints, that further improve the Bayesian mixture model by reducing the proportion of false positive findings among the identified interactions and successfully identifying the reported associations. This is practically appealing for the study of investigating the causal factors from a moderate number of candidate genetic and environmental factors along with a relatively large number of interactions. The natural and orthogonal interaction (NOIA) models of genetic effects have previously been developed to provide an analysis framework, by which the estimates of effects for a quantitative trait are statistically orthogonal regardless of the existence of Hardy-Weinberg Equilibrium (HWE) within loci. Ma et al. (2012) recently developed a NOIA model for the gene-environment interaction studies and have shown the advantages of using the model for detecting the true main effects and interactions, compared with the usual functional model. In this project, we propose a novel Bayesian statistical model that combines the Bayesian hierarchical mixture model with the NOIA statistical model and the usual functional model. The proposed Bayesian NOIA model demonstrates more power at detecting the non-null effects with higher marginal posterior probabilities. Also, we review two Bayesian statistical models (Bayesian empirical shrinkage-type estimator and Bayesian model averaging), which were developed for the gene-environment interaction studies. Inspired by these Bayesian models, we develop two novel statistical methods that are able to handle the related problems such as borrowing data from historical studies. The proposed methods are analogous to the methods for the gene-environment interactions on behalf of the success on balancing the statistical efficiency and bias in a unified model. By extensive simulation studies, we compare the operating characteristics of the proposed models with the existing models including the hierarchical meta-analysis model. The results show that the proposed approaches adaptively borrow the historical data in a data-driven way. These novel models may have a broad range of statistical applications in both of genetic/genomic and clinical studies.

Rapid assessment of genotype-by-environment interactions and heritability for growth rate in aquaculture species using in vitro fertilisation and DNA tagging

Commercial environments may receive only a fraction of expected genetic gains for growth rate as predicted from the selection environment This fraction is the result of undesirable genotype-by-environment interactions (G x E) and measured by the genetic correlation (r(g)) of growth between environments. Rapid estimates of genetic correlation achieved in one generation are notoriously difficult to estimate with precision. A new design is proposed where genetic correlations can be estimated by utilising artificial mating from cryopreserved semen and unfertilised eggs stripped from a single female. We compare a traditional phenotype analysis of growth to a threshold model where only the largest fish are genotyped for sire identification. The threshold model was robust to differences in family mortality differing up to 30%. The design is unique as it negates potential re-ranking of families caused by an interaction between common maternal environmental effects and growing environment. The design is suitable for rapid assessment of G x E over one generation with a true 0.70 genetic correlation yielding standard errors as low as 0.07. Different design scenarios were tested for bias and accuracy with a range of heritability values, number of half-sib families created, number of progeny within each full-sib family, number of fish genotyped, number of fish stocked, differing family survival rates and at various simulated genetic correlation levels

Spoilage pattern of five species of Australian prawns: Deterioration is influenced by environment of capture and mode of storage

Five species of commercial prawns Penaeus plebejus, P. merguiensis, P. semisulcatus/P. esculentus and M. bennettae, were obtained from South-East and North Queensland, chilled soon after capture and then stored either whole or deheaded on ice and ice slurry, until spoilage. Total bacterial counts, total volatile nitrogen, K-values and total demerit scores were assessed at regular intervals. Their shelf lives ranged from 10-17 days on ice and >20 days on ice slurry. Initial bacterial flora on prawns from shallower waters (4-15m) were dominated by Gram-positives and had lag periods around 7 days, whereas prawns from deeper waters (100m) were dominant in Pseudomonas spp. with no lag periods in bacterial growth. The dominant spoiler in ice was mainly Pseudomonas fragi whereas the main spoiler in ice slurry was Shewanella putrefaciens. Bacterial interactions seem to play a major role in the patterns of spoilage in relation to capture environment and pattern of storage

Crop and environmental attributes underpinning genotype by environment interaction in synthetic-derived bread wheat evaluated in Mexico and Australia

Synthetic backcrossed-derived bread wheats (SBWs) from CIMMYT were grown in the north-west of Mexico (CIANO) and sites across Australia during 3 seasons. A different set of lines was evaluated each season, as new materials became available from the CIMMYT crop enhancement program. Previously, we have evaluated both the performance of genotypes across environments and the genotype x environment interaction (G x E). The objective of this study was to interpret the G x E for yield in terms of crop attributes measured at individual sites and to identify the potential environmental drivers of this interaction. Groups of SBWs with consistent yield performance were identified, often comprising closely related lines. However, contrasting performance was also relatively common among sister lines or between a recurrent parent and its SBWs. Early flowering was a common feature among lines with broad adaptation and/or high yield in the northern Australian wheatbelt, while yields in the southern region did not show any association with the maturity type. Lines with high yields in the southern and northern regions had cooler canopies during flowering and early grain filling. Among the SBWs with Australian genetic backgrounds, lines best adapted to CIANO were tall (>100 cm), with a slightly higher ground cover. These lines also displayed a higher concentration of water-soluble carbohydrates in the stem at flowering, which was negatively correlated with stem number per unit area when evaluated in southern Australia (Horsham). Possible reasons for these patterns are discussed. Selection for yield at CIANO did not specifically identify the lines best adapted to northern Australia, although they were not the most poorly adapted either. In addition, groups of lines with specific adaptation to the south would not have been selected by choosing the highest yielding lines at CIANO. These findings suggest that selection at CIMMYT for Australian environments may be improved by either trait based selection or yield data combined with trait information. Flowering date, canopy temperature around flowering, tiller density, and water-soluble carbohydrate concentration in the stem at flowering seem likely candidates.

Classifying genotype by environment interactions for targeted germplasm deployment with a focus on Eucalyptus

A novel methodology for describing genotype by environment interactions estimated from multi-environment field trials is described and an empirical example using an extensive trial network of eucalypts is presented. The network of experiments containing 65 eucalypts was established in 38 replicated field trials across the tropics and subtropics of eastern Australia, with a selection of well-tested species used to provide a more detailed examination of productivity differentials across environmental gradients. By focusing on changes in species’ productivity across environmental gradients, the results are applicable for all species established across the range of environments evaluated in the trial network and simultaneously classify species and environments so that results may be applied across the landscape. The methodology developed was able to explain most (93 %) of the variation in the selected species relative changes in productivity across the various environmental variables examined. Responses were primarily regulated by changes in variables related to water availability and secondarily by temperature related variables. Clustering and ordination can identify groups of species with similar physiological responses to environment and may also guide the parameterisation and calibration of process based models of plant growth. Ordination was particularly useful in the identification of species with distinct environmental response patterns that would be useful as probes for extracting more information from future trials.